Renormalised nonequilibrium quantum field theory: scalar fields 
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We discuss the renormalisation of the initial value problem in quantum field theory using the two- 
particle irreducible (2PI) effective action formalism. The nonequilibrium dynamics is renormalised 
by counterterms determined in equilibrium. We emphasize the importance of the appropriate choice 
of initial conditions and go beyond the Gaussian initial density operator by defining self-consistent 
QQ . initial conditions. We study the corresponding time evolution and present a numerical example 

' which supports the existence of a continuum limit for this type of initial conditions. 
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I. INTRODUCTION 



Nonequilibrium field theory is receiving an increasing level of attention from the side of cosmologists as well as from 
the heavy ion community. The reheating of the postinfiationaxy universe [l| the dynamics of symmetry breaking 
and the formation and decay of cosmological defect networks Q as well as the phase transitions in the early universe 
^ • Q with possible relic gravitational waves [5| are just some of the examples that require the study of out-of-equilibrium 
fields in a cosmological context. Similarly, the rapid thermalisation of the hot quark-gluon plasma 0,0] and its driving 
^ . mechanisms, such as the Weibel instability Q raise questions in the realm of nonequilibrium field theory. 

One of the simplest and most popular strategies to describe an out-of-equilibrium field theory is the classical ap- 
proximation. It has been extensively used for reheating models of the early Universe ^ and also for predicting the 
corresponding production of gravitational waves [l^]- Other cosmological applications include electroweak baryo- 
' genesis [ll|, as well as theories accomodating non- fundamental strings n§\ or domain walls fl3, ^l^l, where the 
nonperturbative treatment is essential. The extreme excitation of the giuon field in a heavy ion collision has also 
made the classical strategy applicable and very successful at early times after a collision [3, |T3| . 
. The success of classical field theory indicates that many of the interesting phenomena in high energy physics are 
actually classical. Indeed, genuine quantum effects play little role if the classical modes are highly excited, and these 
classical fluctuations can play the role of quantum particles. In order for the classical approach to work the UV 
modes must remain unexcited to avoid Rayleigh-Jeans divergences. In fact, this restricts the classical treatment to 
far-from-equilibrium settings, and the classical dynamics automatically drives the system out of its range of validity 
in the course of equilibration. 

It is still possible to split the momentum space of a theory into different momentum regions, where hard degrees 
of freedom follow a quantum-mechanically correct Hard Thermal Loop (HTL) dyiiamics to some finite perturbative 
?H ' order, while the infrared part follows the classical non-perturbative dynamics [1^, [l^ . This allows a kinetic description 
5^ 1 [20I . [2l| for the hard modes in terms of a Vlasov equation . The interplay between classical waves and particles can 
give account for non-trivial dynamics, such as the development of plasma instabilities [H, In this way one can 
avoid the problem of ultraviolet divergences, but the scale separation is not always natural, especially if the coupling 
is not small. 

A step towards the inclusion of quantum corrections from first principles is the Hartree approximation j25l . [26j . It 
assumes a constantly Gaussian density operator and allows to account for quasi-particles propagating in arbitrary 
inhomogeneous backgrounds, which can be as complicated as a network of topological defects [23, Although it 
completely neglects the scattering of the quasi-particles on each other, non-perturbative particle creation mechanisms, 
like tachyonic instability 29] or parametric resonance [30| are within its range of validity. Renormalisation in this 
framework has already been discussed at length [3l[. Despite of its simplicity and clean formulation the fact that 
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thermalisation can not be described in this way [32l . l33l | explains why this approximation scheme could not reach a 
wide acceptance. 

The two-particle irreducible (2PI) effective action provides a first principles approach to quantum field theory 
[3^ [35}. The systematic approximations, obtained from any small parameter expansion of the 2PI effective action, 
usually resum an entire series of ladder diagrams which in turn play a particularly important role in solving the 
secularity problem of out-of-equilibrium perturbation theory. This type of resummation is also present in the Kadanoff- 
Baym equations ^36*1 as well as in Boltzmann equations [s^ . It has been shown that the first nontrivial truncation 
of the 2PI effective action beyond two-loop order already provides a sufficient framework to describe thermalisation 
in scalar theories [H, [H, HO, |4l|, \^ as well as in a model with fermions [i^ . The 2PI effective action has become 
a standard framework for nonequilibrium quantum field theory (4^ . at present, with mostly scalar applications of 
cosmological interest [il,|46, 47]. 

The success of these practical applications has also encouraged more formal investigations on the very foundations 
of the 2PI (and more generally nPI [i^) approach. To be considered as a sensible approach, the latter should 
lead to approximations which reflect, as much as possible, the basic properties of quantum field theories. The 
thermodynamical consistency as well as the energy conserving nature of the time evolution have already been known 
|49l|. U nderstanding how global and local symmetries appear at the level of the 2PI effective action has been studied 
in [set [Slj . Because some important aspects of the dynamics of hot gauge t heories require a 3PI analysis [13, , 
applications have also shifted focus to higher nPl effective actions [s^ ISSIISGI. [57| . For low orders, however, (e.g. for 
the truncation used in this paper) any higher nPI effective action yields the equation of motion what one also finds 
in 2PI 0. 

Another imp ortant issue is renormalisation. So far, the latter has been considered in equilibrium for scalar (sol . 
[ssl . [sol . [60I [6II [6^ . l63j . fermionic [U as well as abelian gauge fields [g^. However, no similar studies exist so far 
out-of-equilibrium beyond Hartree approximation, and most applications were based on cut-off theories defined in 
terms of bare parameters. In this work, we intend to fill this gap and show how renormalisation results obtained in 
equilibrium could be used out-of-equilibrium. For illustration, we consider a scalar {X/4\)(p'^ theory but our approach 
could be extended to other theories of relevance. 

In Secini we recall basic definitions and results in equilibrium. Special attention is paid to the equivalence between 
different contours in the truncated 2PI framework, which we use later on. In Sec. lIIII we move on to out-of-equilibrium 
situations and point out the relevance of beyond-Gaussian initial conditions for continuum field theories. We suggest a 
self-consistent initial condition and study the corresponding time evolution in the 2PI three-loop approximation. Our 
numerical results strongly suggest the existence of a continuum limit for such initial conditions, unlike what happens 
with Gaussian initial conditions. Appendices [Cl and [P] collect the algorithms used to obtain the results in Sec. IIIII 



II. EQUILIBRIUM QUANTUM FIELD THEORY 

In this section we consider a real scalar field in equilibrium at a temperature T = 1//3. In this context, one is usually 
interested in determining thermal expectation values of products of field operators ordered along a given time contour 
C. In Sees. Ill AIIII Cl we recall the definition of the contour-ordered propagator in equilibrium and explain how to 
evaluate it non-perturbatively within the 2PI approximation scheme. This requires non-perturbative renormalisation 
as we discuss in Sec. IIIDI 



A. Contour-ordered propagator 

Consider a complex time contour C, the properties of which we shall specify in what follows. The contour-ordered 
propagator is defined as 

G{x,y) = Qc{xo,yo)G>{x,y) + ec{yo,xa)G>{y,x), (I) 

where Qci^o, yo) denotes the step function along the contour, equal to one if yo precedes xo and zero in the opposite 
case. The Wightman function[70j G^{x,y) appearing in Eq. ([T|) is defined as 

G^{x,y) = {ip{x)ip(y))p = ^f~m • (2) 

In equilibrium G^{x,y) depends on the difference of its arguments only and we shall use this simplification 
when necessary. Using Lehmann's representation one shows that {x — y) is analytic in the complex domain 
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—(5 < Im (xq — hq) < 0. From the equal-time commutation relations, we have 

[G> [x, y) - G> (y, x)] = ~i6^''^ [x - y) . 
Finally, one can easily show the KMS (Kubo-Martin-Schwinger) relation 

G> {xo - i/3, yo- x, y) = G> (yo, xq] y, x) , 
which holds for < Im (sq — yo) < In Fourier space, [7l| it reads 

G>{-po,-p) = e-^P"G>ipo,p). 



(3) 
(4) 

(5) 
(6) 



From the analyticity property of G^ {x, y) it follows that the contour-ordered propagator G(x, y) is properly defined 
only for contours with a decreasing imaginary part. Moreover, when one tries to practically compute G(x, y), a second 
important constraint appears on the contour, namely that it should stretch from some - arbitrary - initial time ti to 
a final time tj — i/3. This condition, which appears both in the canonical and path integral approaches, is dictated by 
the presence of the operator e~^^ in the thermal average of Eq. The same condition has been recently found for 
the convergence of real-time lattice simulations (66j . One can construct numerous contours with decreasing imaginary 
part and stretching from time tj to tj — i(3. We shall refer to them as admissible contours. If one is interested in 
real-time aspects, the first branch of the contour should be contained in the real time axis. A particular example of 
such a contour is the close-time path, represented in Fig. [T] 



FIG. 1: Close-time path. We here chose ti = 0. 



B. Propagator from the 2PI approximation scheme 

Any propagator ordered along an admissible contour admits the path integral representation 

J Vip ip{x)ip{y) exp {ij^£{x)} 



G(x,y) 



JVip exp {i J^£{x)} 



where £{x) is the Lagrangian density 



C{x) = lidpix)f ~ \m'i^ix)f - ^Mx)r 



(7) 



(8) 



and the path integral is taken over fields satisfying the boundary condition ^{ti; x) = ip{ti — i(3; x). This path integral 
representation can be used to obtain a perturbative expansion of the contour-ordered propagator. It can also be used 
to define systematic non-perturbative approximations of the latter, as we now explain. 

The contour-ordered propagator ([7]) can in fact be obtained from the generating functional 



W[J,K] = -i\n / Vif e'Jc 



(9) 



either from a second derivative with respect to the source J(x) or a single derivative with respect to the source K{x, y), 
evaluated for vanishing sources. If we now introduce, the conjugate variables (t>{x) and G{x, y) defined by 



6W[J, K] 
SJix) 



(x) and 



SW[J, K] 
SK{x,y) 



G(x,y) + cf>{x)c^{y) 



(10) 
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we can construct the double Legendre transform of W[J, K]: 



, G] ^ W[J, K] - I J{x)c^{x) -\n K{x, y) G{x, y) + ^{x)^{y) , (11) 



'c JC ^ 

from which one can easily recover the initial sources J(x) and K(x, y) as 



- -J{x)- K{x,y)(t){y) and —. = --K{x,y) . (12) 



Scj^ix) ' ' Jc 5G{x,y) 2 

The benefit of this last equation is that, for a given source K{x,y), it defines the two-point [t^ function G{x,y) as 
the solution to an implicit variational equation, which in turn is a powerful means for resumming Feynman diagrams 
and defining the contour-ordered propagator non-perturbatively. 
It is particularly convenient to introduce the decomposition 

T[^,G]=So[^]+'-Tyc lnG-i + ^TrcGoiG + ri„t[0,G] (13) 

where So[ip] denotes the free classical action along the contour C and G^^(x,?/) = i{d^ + m'^)Sc{x,y) is the corre- 
sponding free inverse propagator. This decomposition allows one to rewrite Eq. (|12p as 

5c{x, y)= I d^z [G^\x, z) - S(x, z) - iK{x, z)] G{z, y) , (14) 
Jc 

where the self-energy S(a;, y) is obtained from 

5G{y,x) ■ 

Since iT-mt [4>i G\ has a simple diagrammatic interpretation [35l | as the sum of all vac-to-vac two-particle-irreducible 
(2PI) diagrams of the shifted interaction part of the theory Sirit[4' + (^s], it is very easy to define systematic non- 
perturbative approximations of the contour-ordered propagator using Eqs. (jl4p and psp . One has simply to select a 
certain number of 2PI diagrams in rint[0, G], plug them in the right-hand-side of Eq. (fTS)) and solve the latter together 
with Eq. (fT4|) . In this work we concentrate on the three- loop approximation to Finti^, G] in the absence of a field 
expectation value for which Eq. becomes 

S(a;,y) = -i'i^o{x)6c{x,y) + Qc{xo,yo)T:> {x,y) + Qc{yo,xo)T,>{y,x) (16) 

with 

So(a:) = ^G(a;,a;) and I]> (x, y) = G> (x, y)^ . (17) 
2 o 



C. Equivalence of contours 



Because we are discussing equilibrium, we set the source K{x,y) to zero for the moment. We shall later consider, 
in Sec. Illlj a non- vanishing source has a means to bring the system out-of-equilibrium. 

Inverting equation Eq. (fTi|) to obtain the propagator G{x,y) in terms of the self-energy T,{x,y) is usually rendered 
cumbersome due to the non-real parts of the contour. As an important simplification, in App. [X] we prove that, in 
equilibrium, any propagator ordered along an admissible contour can be reconstructed from the so-called real-time 
propagator, obtained from solving Eqs. and on the real-time path (depicted in Fig. (2), with the KMS and 
equal-time commutation relations as boundary conditions. The KMS condition is necessary because, in going from 

FIG. 2: Real-time path. 

the close-time path to the real-time path, one looses any reference to the temperature and one needs a sensible way 
to reintroduce it. 
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In what follows, we restrict our analysis to the real-time path. In this case, Eq. ((T4)) can be put in a more tractable 
form by introducing the functions 



as well as 



F{x,v) 



1 



1 r 



G> (x, v) + G> (y, x) , v)=tG>{x,v)- G> (y, x) 



S^(a;,y) + I]>(y,x) , T.p{x,y) = i {x , y) - J:> {y , x) 



(18) 



(19) 



Using the equal-time commutation relations, one can then recast Eq. (jl4p into a pair of partial differential equations 



{dl + + Eo) F{x) = / dz^^F{x - z)p{z) - / dz'^T.p{x - z)F{z) , 



(20) 



[dl + + Eo) p{x) = - / dz^Y.p{x - z)p(z) , 



(21) 



where Eq, Ei?(x) and Ep(a;) can also be expressed in terms of F{x, y) and p{x, y). In writing these equations, we have 
used the fact that, in the absence of sources the functions F{x,y), p{x,y), Y^pix^y) and Yip{x^y) only depend on the 
difference of their arguments. 

In order to yield a particular solution, the previous two equations need to be supplemented by some boundary 
conditions. Those for the spectral density p{x) are fixed by the equal time commutation relations: p{x)\^^^q = and 
dxoP{x)\^^^Q — 1. As for F{x), the boundary condition is nothing but the KMS condition we alluded to before. In 
terms of F / p components, it takes the form 



(22) 



where /(po) — ^/{e^^° — 1) denotes the Bose-Einstein factor. In fact once the KMS condition is assumed, the system 
of equations ((20l) - pT|) becomes redundant. To sec this in our particular example, notice that 



E>(p): 



A2 



6 J (27r)4 J (27r)4 



dH 



G> {k)G> [l)G> {p - k - I) . 



(23) 



from which one can check that E^(p), and in turn T,p{p) and Ep(p), obey KMS conditions similar to Eqs. ^ and 
(|22| respectively. Then writing Eqs. (pO| and (|2T|) in Fourier space (we set iJp= -j? + rn? -I- Eg): 



dhJ 



{-pi+uji)F{p)= / — 
{^pl+uj})p{p)= / ^ 



T.f{p) p{uj\p) , T,p{uj;p} F{p) 



i{pQ ~ LU - ie) i{po—LLi + ie) 

^p{p)p{'^;p) , ^p{^;p)p{p) 



i{po — uj - ie) i{po—LO + ie) 



(24) 
(25) 



it is straightforward to check that Eq. implies Eq. ([^5]) and vice versa. We conclude that in equilibrium, one only 
needs to solve Eq. (pi]) for p{x) and determine F{x) from the KMS condition. Equation (PT|) needs to be renormalised, 
as we now explain. 



D. Renormalisation 



Equation (^J) for the spectral density is usually plagued by ultraviolet divergences which one needs to renormalise 
in order to define a continuum limit. It has to be noted that, due to the presence of a Landau pole, this limit does 
not exist in the strict sense. It is true, however, that there is a wide range of couplings where the Landau pole is 
far in the UV and renormalisation can ensure the insensitivity to the cut-off, if it is significantly higher than other 
physical scales, but does not exceed the Landau pole, at which the computed bare coupling diverges. We use the term 
"continuum limit" in this restricted sense. Pattern of the divergence of bare parameters in the vicinity of the Landau 
pole has been discussed in f67| in the 2PI three-loop approximation, that we also use here. 
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Renormalisation on the real-time path has been considered in [5^ . In this work we shall rather consider renormal- 
isation in the so-called imaginary-time path and infer renormalisation on the real-time path. Indeed, according 
to App. lAl the solution to Eqs. (fT4|) and (fT6|) on the real-time path supplemented by the KMS and equal-time com- 
mutation relations can be used to construct propagators ordered along any admissible contour, at the same level of 
approximation. In particular, if one chooses the imaginary-time path depicted in Fig. [3l one obtains the so-called 
imaginary-time or Euclidean propagator (— /? < t < /?) 

GE{T;x) = eiT)G>{-iT;x) + e{-T)G>{iT;-x), (26) 

where {—'iT; x) is the analytic continuation of the real-time Wightman function to the imaginary-time path, as 
defined in Eq. (IA3P of App. lAl The KMS condition implies that Ge{t; x) is /3-periodic in r from which one concludes 



t 



FIG. 3: Imaginary time path. 

that Ge(t] x) can equally be represented by its Fourier modes Ge{p) with p = (iujn]p) and tj„ — {2tt / (3)n a discrete 
Matsubara frequency: 



GE{ii^n;p) = J drj d^xe''^"^e-'P-^G>{-iT;x). 

Using Eq. ()A3p and the KMS condition, one then shows that 

d'^p ip{pa;p) 



(27r)4 itUn - Po 



(27) 



(28) 



from which it follows that 



p(Po;p) = 2iImG£;(po + • (29) 

We conclude from this, that in order to renormalise p{pQ;p), it is enough to renormalise GE{i^n]P)- This assumes 
that the analytic continuation involved in Eq. (j29D does not bring in new divergences. This assumption is plausible 
since, in what follows, we restrict to a spatial cut-off and leave the time direction in the continuum, see below. 



1. Imaginary time counterterms 

In Fourier space, the Euclidean propagator in the approximation at hand is such that G^^{p) = G^^g(p) + S£;(p) 
with G^ g(p) = LUj-^ +]f + rn? and 

I]b(p) = 5Y.%(p) + \fGE{k) - y GE{k)GE{l)GE{p - k - I) , (30) 

The piece 5E^(p) contains counterterms which should be adjusted in such a way that UV divergences are absorbed. 
Due to Eq. (j28p . the counterterms in imaginary-time are exactly those needed to renormalise the equations in real- 
time. Strictly speaking this is true if the time integral is not discretized. For this reason we introduce a spatial cut-off 
only. For such an anisotropic regularisation and at three-loop order in the 2PI-loop expansion the counterterm piece 
reads 

(5I]^*(p) = 6Zt Lul + SZs f + 5M^ (31) 

with 

6M^ = 6ml + ^ j:GE{k). (32) 
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The three counterterms in Eq. (j3ip can be fixed using renormahsation conditions at a reference temperature T* and 
a reference momentum p* 



d^%{p) 



dpi 



= 0, (33) 

p=p* 



see App. [Clfor Algorithm ([T]) which finds these counterterms. These three renormahsation conditions are, however, 
not sufficient to ensure insensitivity to the UV cut-off at a temperature T away from the renormahsation temperature 
T*. For this to be achieved one needs to adjust the couphng counterterm in such a way that couphng UV 
subdivergences which appear in Tie{p) are properly absorbed. These subdivergences can easily be accounted for by 
remarking that they are obtained after opening a perturbative line (corresponding to a free propagator Ge.o) in any 
of the diagrams contributing to Tje{p)- Algebraically, this corresponds to considering the function 5Yie{p) / Efi{<l) 
which is shown to be given by: 

5^e{p) _ VE{p,q)/2 ^^^^ 



SGeAq) (l + G£,o(g)SB(q))2 
with 

VE{p,q)=AE{p,q) -^j:VE{p,r)[GEir)fAEir,q) (35) 

and 

Ae{p, q) = X + SXE- \^jGE{k)GE{p -q-k). (36) 

Equation (j34p tells that if one wants to properly absorb coupling subdivergences appearing in T,E{p)^ one needs to 
renormalise the function Ve (p, q) as well. This is actually done by means of the following renormahsation condition 
at the reference temperature T*: 

V^{p\pn^\. (37) 

This fixes the value of the coupling counterterm SX e and ensures that Ve {p, q) as well as (p) are not sensitive to 
the UV cut-off at any temperature T. 

Although it is possible to solve Eq. (I35p exactly as a linear algebra problem, we recommend the strategy detailed 
in Algorithm see App. [Cl In Ref. 67] we have solved this equation numerically in imaginary time and calculated 
the renormalised 2PI pressure of a scalar (p"^ theory to three-loop order in the 2PI-loop expansion. 



2. Real-time counterterms 



In principle the counterterms in real time should be equal to those in imaginary time, as long as time is not 
discretised. We have checked this numerically for (p^ theory, but the agreement was not accurate, unless we used 
extremely anisotropic regulators {at /a < 0.05), which is a highly inconvenient choice for the subsequent nonequilibrium 
application due to the expense of the required storage and computation. If instead, we want to use a conveniently 
coarse discretisation, we should not rely on the counterterm values obtained in imaginary time but rather compute 
them in real time, directly. Still, calculations in imaginary time are not totally useless here since they allow for a 
cheap determination of real time counterterms (without actually solving a real time version of Eq. (|35p ). as we now 
explain. 

First of all, we notice that in imaginary time the field strength counter-terms 5Zt and 5Zs grow logarithmically 
with the cutoff with a prefactor of the order of lO^'^. This means that the finite part of the self-energy reaches a 
cut-off insensitive value before field strength divergences can even show up. In other words, for the cut-off values we 
use, we can equally drop 5Zs and 5Zt and still obtain cut-off insensitive results. [tII We assume that this is also true 
in real time and check this a posteriori. 

We still need to fix the mass and coupling counterterms in real time. The trick we use here to avoid solving the 
real-time version of Eq. (j35p is that the coupling can equally be fixed by conside ring any other coupling dependent 



observable, such as the thermal mass or the renormahsation scale dependence [61[. Suppose then that we know 
the thermal mass mjj^(r) from a calculation in imaginary time (we may use here a convenient close-to-continuum 
regulator, e.g. at <C a, aT <C 1, am <^ 1), defined as 

mt\(r) = S£(p = 0), (38) 
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where Tie{p = 0) contains the imaginary time counterterms Sm^ and SXe, and renormaHsation is implemented accord- 
ing to Eqs. and ([57)) . see Algorithms ^ and We can evaluate this thermal mass at different temperatures. 
In particular if we evaluate it for two different temperatures we have in principle enough information to recover the 
values of the counterterms Srn^ and 5Xe- Now, the thermal mass has also an expression in real-time which involves 
the real-time counterterms 

= Srn' + j J^^Pi^) + dtJ:,{t;p = 6) . (39) 

Using the two values for the thermal mass obtained in imaginary-time we can obtain the values for the real-time 
counterterms. This is used in Algorithm (3), see App. [Cl 



III. DEPARTING FROM EQUILIBRIUM 



Let us now consider a nonequilibrium situation. If we aim at studying late-time dynamics and thermalisation, the 
evolution should include the counterterms obtained in the previous section in order to describe the final equilibrated 
state properly. Because these counterterms do not depend on time, they will be present at any time of the evolution. 
Then, depending on how one sets up the evolution equations, this might lead to problems such as a UV divergent time 
evolution. In Sec. 3.1, we first quickly revisit the evolution equations with a Gaussian initial condition and discuss 
why one fails in obtaining a continuum limit for this type of evolution. In Sec. 3.2, we start the evolution from a 
self-consistently dressed initial quantum state. We give arguments for the existence of a continuum limit in this case 
and support them with numerical evidence. 

A. Gaussian initial condition 

A popular initialisation of the propagator is a Gaussian quantum state (39j . A remarkable feature of this special 
choice is that the memory integrals start at initial time, and the only input parameter to be fixed is the initial 
propagator and its first derivatives at equal time. The equations of motion in the F/p formalism read [39|: 

ya Xo 

{dl + m^+j:oix))F{x,y)=Jd'^zJ:^{x,z)p{z,y)-Jd*z^P{x,z)F{z,y), (40) 



Xo 

{dl + m^+j:o{x)) p{x, y)^ - j d^zJ:P{x, z)p{z, y) . (41) 

yo 

This initial prescription has been frequently used to successfully describe systems with a fixed cut-off. In low 
dimensional systems renormalisation is less problematic [6^ . In a 3-1-1 dimensional continuum quantum field theory, 
however, this class of initial conditions is not useful for the continuum limit of the time evolution is not accessible. 
One simple way to view this is to consider the equation for F(x^ y) at initial time xq = yo = 0. The contribution So (2;) 
in the left-hand-side of this equation contains the counterterms which are supposed to renormalise the equilibrium 
state obtained from the time evolution. In contrast, the right-hand-side of the equation is equal to zero and thus, 
part of the diagrams which should be absorbed by the counterterms are absent: One has an unbalanced divergence. 

In Fig. m we present numerical results for the time evolution obtained from this kind of initial condition. We 
parametrise the initial equal-time statistical correlator F{t,t;p) using the particle number 

nip) = NeM-m-Pcf/^a^) (42) 

with N = 5, a = 0.6m, Pc = m. We use the counterterms determined at T* — m and plot the evolution of three 
different modes \p\ = 0.4,0.8, 1.6 for three different values of lattice spacing am = 1/4, 1/6, 1/8. On such fine lattices 
one expects at least an approximate convergence as am — > 0, as it will indeed happen in the next subsection. These 
curves, however, show no sign of a continuum limit. In what follows, we propose a possible way to cure this. 



B. Self-consistent initial conditions 



The general framework of Sec. [TTl suggests an alternative way to bring the system out of thermal equilibrium. One 
simply has to calculate the 2PI propagator in presence of a non- vanishing external source K{x,y) living on the real- 



9 



1000 



100 



10 



0.1 




0.01 



0.1 



10 



100 



time [m"^] 



FIG. 4: Time evolution of the equal-time statistical correlator F{t,t;p) with fixed Gaussian initial propagator (A = 24, 
box size : L = 32) for three different modes \p\ — 0.4, 0.8, f .6 and three different values of lattice spacing am — 1/4, f/6, f/8. 
The counterterms we determine at temperature T* — m are not enough to ensure a continuum limit of the time evolution. 



time path. In this paper we assume the components Kq{x) and Kp{x, y) to vanish, in such a way that the equations 
for the 2PI propagator read 



ya Xo 

{dl+m^ + Y.^{x))F{x,y)J dz^Y.^{x, z)p{z,y)-( dz'^T.p{x, z)F{z,y) 



Xo 

{dl+m'^ + T,Q{x))p{x, y) = - dz'^T,p{x, z)p{z, y) , 



(43) 
(44) 



where Sf^ = ^f +iKp. The benefit of this approach is that the renormalisability of the evolution equations only 
depends on the properties of the source K(x,y). As we motivate in App. |B1 once the problem has been properly 
renormalised in the absence of source, see Sec. IIIDl introducing a source does not bring new divergences provided 
the UV behavior of the source is well under control. More precisely we show that, on the imaginary-time path, a 
source K with proper asymptotics does not alter the UV structure of the Euclidean theory. To do so we expand the 
self-energy E^; in powers of K around the equilibrium solution which we know how to renormalise: |74| 



AS, 



z — / „ 



n>l 



! 



(45) 



K=0 



and show that each term of this expansion is finite provided K has the correct UV asymptotics. Of course, strictly 
speaking, it is not completely obvious that the result of App. [B] obtained on the imaginary-time path can be applied 
to the real-time path. We shall however provide numerical evidence for a continuum limit on the real-time path which 
is an indication that our argument does not really depend on the contour we consider. 
Suppose we use a source given by 



KAx,y)^l^^^''-y^ if2;„<0 and 
^ ' [0 if > or 



yo < 0; 
or yo > 0. 



(46) 



For times smaller than zero, the system is in a steady state sustained by the translationally invariant source Kp{x — y). 
Such states exist for all times and can be obtained by solving the simplified, translationally invariant equations 



u 

]o) Fix)=Jc 



= I dz^Sf (x - z)p(z) - I dzJj:p{x - z)F{z) 



(47) 
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{dl + + So) pix) = - j dz^Y.p{x - z)p{z) . (48) 



The source Kp[x — y) alone does not completely fix the solution of these equations, as it was also the case in 
equilibrium (no source). To define a particular solution we introduce a non-thermal boundary condition in the form 
of a generalised KMS condition 

F{p^-p) = Q + /(Po;p)) p(Po;p) , (49) 

where f{poiP) is an arbitrary function. The ad-hoc use of Eq. has already been suggested in Ref. fiS] as a recipe 
to prepare a dressed initial state. We shall restrict to this type of boundary conditions in what follows. Notice that 
this generalised KMS condition actually fixes the source to 

Kf{p) = - + /(Po;p)) Sp(p) + I ^f{p) , (50) 

as it can be shown by writing Eqs. (|47p - (l48p in Fourier space. We need to require that f{po\p) approaches a thermal 
Bose-Einstein factor with some temperature sufficiently fast in the UV. This restriction makes sure that the expansion 
in K used in the above argument is defined around equilibrium and that the asymptotics of K is appropriate to ensure 
that the steady solution admits a continuum limit. The details on how to obtain the steady propagator are given in 
Algorithm (4), see App.O 

Now let us turn to the nonequilibrium dynamics. For times greater than zero, one has to solve 

ya Xo 

{dl+m^ + ^o)F{x,y)^ ( dz^^F{x,z)p{z,y) - ( dz^^p{x, z)F{z,y) (51) 



{dl + w? + llQ)p{x,y) = - jdz^T.p{x,z)p{z,y) , (52) 



ya 

where in the integrals, the parts involving times smaller than zero involve the pre-calculated steady solution we just 
discussed. In the course of time evolution xq is always considered as the most recent time. Notice, that for practical 
purposes, we only keep the latest part of the memory integrals with jxo — ^ol < imom- Equations ([5T|) - ([5^ look similar 
to (|40p and (|¥T|) but there is in fact a big difference: The memory integrals are present at time xq = yo — Q already 
and can therefore prevent an unbalanced divergence, unlike what happened in the case of a Gaussian initial condition. 
Still for the divergences to be completely cancelled by the counterterms, the inhomogeneous source K, see Eq. (j46l) . 
should again have the required asymptotics. This in turn depends on how one chooses f{pQ;p) in the generalized 
KMS condition. In this paper we shall not prove analytically that this source has the correct asymptotic behavior. 
Rather we shall provide numerical evidence for this, by generating a cut-off insensitive time evolution. 

In Fig. [5] we present the time evolution of the equal-time statistical propagator F{t,t;p) for three different modes 
\p\ = 0.4, 0.8, 1.6 and three different values of lattice spacing am = 1/4, 1/6, 1/8. We prepare the initial correlated 
state with 

where T{p) — T* + cjp>/log(l -I- l/n{p)), uj^ = + m^^^{T*) and n{p) is the initial particle distribution (|^ that we 
used in the case of the Gaussian initial condition. The counterterms are again determined in equilibrium at temper- 
ature T* — TO. For a given mode, the curves representing runs for different values of the spatial cut-off lie almost 
exactly on top of each other, showing that the continuum limit has been reached. This plot strongly suggests that 
our approach to out-of-equilibrium renormalisation correctly removes all UV divergences and thus allows to define a 
continuumlimit . 

The final relaxation to equilibrium is particularly sensitive to the truncation of the memory integral. For the plot 
we used tmcm = 12to^^. Keeping longer memory we could avoid the small deviation of the curves at ~ lOOOm"^. 
Long memory integrals, however, are vulnerable to "instabilities" on the course of the time evolution. Traces of such 
"instabilities" can already be seen in Fig. [5] on the uppermost curves at ~ IOOto^^. This phenomenon manifests 
in a high-frequency oscillation on the low-momentum F{t,t;p) curves, and eventually relaxes. Also, after this high- 
frequency oscillations have died out, the F(t,t;p) function continues normally and maintains the approximate am- 
independence. This phenomenon persisted in an increased volume. Nevertheless, from our test run with L = 48 we 
can conclude that the curves in Fig. [5] are very close to the infinite volume limit up to about ^ IOOOto"^ where finite 
volume effects start to matter. 
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FIG. 5: Time evolution of F(t,t;p) (A = 24, box size : L — 32) for three different modes \p\ — 0.4,0.8, 1.6 and three different 
values of lattice spacing am = 1/4, 1/6, 1/8. The counterterms are determined at temperature T* — m. For a given mode, 
the curves corresponding to runs with a different lattice spacing almost exactly match each other, which indicates that the 
continuum limit has been reached. We have checked that the curve for am = 1/2 still shows sizeable deviations. 



We have discussed the renormalisation of the initial value problem in a scalar quantum field theory. We have 
introduced a self-consistent initialisation of the 2PI equations of motion that defines a cut-off independent initial 
quantum state. We have also argued for the absence of new divergences in this setting and performed a numerical 
analysis to show the continuum limit evolution of a nonequilibrium quantum field theory. 

In this paper we considered the excitation of an originally equilibrium system using a two-point source. Although 
we have not discussed it in this paper, one could give similar arguments for the finiteness of the solution in presence 
of a one-point source J(x) with proper asymptotics. For a spatially homogeneous initial condition, J{x) may only 
depend on time. This kind of initialisation is for example suited to describe the parametric resonance scenario. 

The primary goal of this paper is to provide recipes. Although scalar fields have already been extensively discussed 
out-of-equilibrium in the 2PI approach, most algorithms and technical details have not been well documented yet. To 
encourage research groups to start working in this field we provide the detailed algorithms used in this work. 

Although scalar theories are interesting on their own (e.g. in cosmological scenarios), they are often considered 
as a warm up exercise for the more complicated gauge theories. The renormalisation of the 2PI effective action for 
quantum electrodynamics has been understood by now in equilibrium [6^, and a similar nonequilibrium discussion 
and numerical analysis is possible. One issue of interest which could be then studied is the gauge-fixing dependence 
of the time evolution. 
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APPENDIX A: EQUIVALENCE OF CONTOURS 



We would like to show here how to construct any contour-ordered propagator from the real-time propagator defined 
as the solution G{x, y) to Eqs. (jl4p and (jl6p on the real-time path, obeying the equal time commutation relations 
and the KMS condition. As we have seen in Sec. Ill C[ these boundary conditions specify a unique solution. It will 
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be convenient to consider the corresponding Wightman function G^{x, y), defined, for the moment, for real values of 
the time arguments. Combining Eqs. and ([^5]) . one obtains an equation for G^{x,y) in Fourier space: 

(-pg + c.|)G>(p) = I ^ 

where ujl = + m'^ + T^o and Sq = (A/2)G(0) = (A/2)G>(0). Using the KMS condition in the form 
p{ijj]p) = ig{uj)G^{ijj;p), where g{uj) = 1 — e^^'^ , and the fact that a similar condition holds between the self-energies 
'Sp{Lu;p) and Y,'^{uj;p), see Sec. Ill Cl we arrive at 



i{po -u) -ie) 



^p{LO-p)G>{p) 



iiPo 



(Al) 



S>(p)G>(c.;p) ^ i:>{u-p)G>{p) 
po — uj — ie po — uj + ie 



(A2) 



We can now use G^{p) to define G^{x, y) for complex values of the time arguments: 

G>{x,y) ^ I ^e-^(--^)G>(p). (A3) 

If we assume that G^ {p) does not grow exponentially at large positive frequencies, this definition makes sense provided 
that — /3 < Im (a;o — yo) < 0. Indeed, if Im [xq — yo) < the integral is well defined as po —> +00. Moreover, thanks 
to the KMS condition which relates the behavior at large positive and negative frequencies, it is easy to see that 
the condition —(3 < Im(a;o — yo) guarantees that the integral is well defined as po —>■ —00. Thanks to the analytic 
continuation (jA3p . we can now define, on any admissible contour C, the following function 

Gc{x,y) = ec{xo,yo)G>{x,y) + ec{yo,xo)G>{y,x). (A4) 

The purpose of such a definition will become clear in a moment. Notice that, due to the equal-time commutation 
relations, acting twice with dx on Gc{x,y) generates a delta function on the contour. Keeping this in mind, one 
arrives at (we consider here the case Iuixq < Im?/o) 

t{dl+cul)Gc{xo,yo;p)=Sc{xo,yo)+ J ^e~^P'>^^'^-y'>\{~pl + u;l)G> (p) . (A5) 

Using Eq. (|A2[) and a simple change of variables, one finally obtains 
Sc{xo,yo) =i{dl„ +1^^) Gc{xo,yo;p) 



g-jc^i(xo-yo)g(^2) - e-*'^2(^o-2'")g(wi) 

X -, ^ , (A6) 

- UJ2) 

where we have dropped the e in the denominator since the integrand is well defined as coi and uj2 become equal. 

Now we would like to show that Gc{xo,yQ;p) is in fact the contour ordered propagator in the approximation at 
hand. In other words, we need to show that Gc fulfills Eqs. (fH)) and on the contour C. To this purpose let us 
evaluate the right-hand-side of Eq. (fTi]) . It can be written as 



R = 'i-{dx„+^p)Gc{xo,yo;p)+ I dzoJ:>{xo,zo;p)G>{yo,zo;-p) 

dzo S > (xo , zo ; G> (zo , 2/0 ; p) 



yo 

ti-iP 

dzaT,^{zo,XQ;-p)G>{zo,yo;p) , (A7) 

2:0 

up — ,,u I So with So = (A/2)Gc(0) and S>(x, y) = — (A^/6)G>(a;, yY . Notice that So coincides with the one 
in Eq. (55) because Gc{x,x) — G{x,x). Moreover, because S^(a;,y) is simply the cube of G-^{x,y), one can use the 
fact that, for any complex times xq and j/q such that — /? < Im (xq ~ yo) < 

S>(x,2/) ^ I ^e-^'(--^)S>(p), (A8) 
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to arrive at 

R = ^i9lo +^f)Gc{xa,ya;p) 



i{uji - UJ2) 

- LO2) 

i(uJi-W2)(ti-il3) _ i(uJx-uJ2)^o 

Y>{-ui: -p) G>{uj2]p)e-'^'''°e''^^y° ^ . 

(A9) 

Using the KMS condition in the form (— Po ; —p) = e~^^°G^ {po ; , it is possible to massage the following expression 
and obtain the right-hand-side of Eq. (|A6p . Thus R = Sc{xo,yo) which proves that Gc is actually the propagator 
ordered along the contour C, as announced. 

APPENDIX B: RENORMALISATION WITH A SOURCE 

Here we would like to study the influence of a source K on the renormalisation of the 2PI self-consistent equations 
for the propagator. To simplify the discussion we consider here the 2PI effective action on the imaginary-time path. 
In presence of a source, the Euclidean propagator is given by (space variables implicit) 



" dr'Gsir, t') iG^^oi^', r") + Sis(r', r") + K{t' , t")] = 5{t - r") (Bl) 



with 

2Sr{ 



int 



5GB(r,r') 



(B2) 



where Ge.o is the free Euclidean propagator and = zFint is the interaction part of the Euclidean 2PI effective 
action. We would like to determine under which conditions the counterterms obtained for K = are enough to 
renormalise in presence of the source K . In order to compare to the situation in the absence of source, we consider 

AI]b(t, r') EE ^E{r, t') - Ef =0(r, r') . (B3) 

Assuming that the amplitude of K is small enough, it makes sense to expand AS ^ as a functional Taylor expansion 
in powers of K . We then would like to find out under which conditions, the different terms of these expansion are 
void of UV divergences. Let us in particular consider the contribution linear in K. At leading order in K , 

AJ^Eir, r') = -l [ VE{r, r; rj, v')G^-°{v, p)G^^°{p\ v)K{p, p') (B4) 

Jr),n',p,p' 



with 



Ve{t,t';X,X')=Ae{t,t';X,X') 
1 



2 



VeIt, t'; ij, 77')Gr°('7, p)Gf =°(p', v')^e{p, p'; A, A') , 



(B5) 



and 



In the previous equations, G^^^, Ae and Ve are evaluated at K = Q, in which case one can use translation 
invariance. It is then convenient to switch to Fourier space and define 

Gf=0(p,p') ^ {27rrS^^Hp + p')GE{p) (B7) 
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and 



AEip,p'; q, q') ^ {2TTf5^^\p + p' + q + q') Ke{p.p' , q) , (B8) 
Ve{p,p'-. g, q') = {2^)H^^\p + p' + q + q') Ve{p,p\ q) , (B9) 



in terms of which, to hncar order in K, 



with 



Ai:E{p,p') = -^t VE{p,p\r)GE{~r)GE{p + p' + r)K{-r,p + p' + r) (BIO) 



VE{p,p',q) = AE{p,p',q) 

- 2 f VEip,p\r)GEi-r)GE{p + p' + r)AE{-r,p + p + r,q) . 



(Bll) 



This last equation is more general than Eq. (|35p. One can, however, show that the counterterm SXe needed to 
renormalise Eq. ([35]) is also the one needed to renormalise Eq. (jBllI) . 75] So no new counterterms are needed here 
and Ve{p,p' ,q) is automatically UV convergent. Since Ve{p,p' , q) grows at most logarithmically with increasing 
momentum q (Weinberg's theorem) and since GE{q) decreases as 1/q^ (up to logarithms), we conclude from Eq. (jBlOp 
that, to linear order in K, ASg is UV convergent provided that, for any value of the fixed momentum u, the source 
K decreases at least like|76| 

Kip,u~p)r.- (B12) 
P 

as p goes to infinity (up to logarithms). A similar analysis shows that, with such a source, higher order contributions 
to AE^; (in powers of K) are also UV convergent. From this we conclude that provided that the source K follows - 
at least - the asymptotic behavior (|B12p . renormalisation in the absence of K is enough to renormalise the system in 
the presence of K. 



APPENDIX C: ALGORITHMS 



The recipes that we disclose here are improved versions of those we used in Ref. [63]. We do not repeat all the 
numerical details that we described there. 

Algorithm 1 Finding the imaginary-time counterterms SZs.t and SM'^ by computing the imaginary-time propagator 
at temperature T* : 

1 . Start from = Gq . 

2. Compute [IttI] the right -hand- side of Eq. JSOD at T = T* keeping the counterterms SZt,s and SM'^ 
as parauneters . 

3. Adjust SZt.s and SAI'^ such that Eqs . 0330 are satisfied and use them to compute an updated 
version G*E^ncw of G*e ■ 

4. Update G^ using G*e = OiG*E ^^^^ + {I — a)G*E where a denotes 
a properly adjusted convergence paramieter between eind 1. 

5. Iterate from step [2] until G^ converges. 

Algorithm 2 Finding the imaginary-time counterterm SXe by solving the Bethe-Salpeter equation for V* (p* . q) -'fT^l 

1. Once and for all, calculate A£;(p, g) at T = T* . 

2. Start from V|(p*,g) = A. 

3. Compute the right hsind side of Eq. g35D for T — T* and p ^ p* , keeping the counterterm 6Xe 
dependence as a parameter. 

4. Solve the linear algebraic equation 0370 for SXe and use the latter to define an uptaded 

VUp\q). 
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5. Iterate from step [3] until V^{p*,q) converges. 

After fixing 6Xe it is not necessary to recalculate with the new counterterms. Neither SZt^s nor SM"^ are 
changed. The value of SXe is important to obtain UV convergent results at temperatures away from T*. 

Algorithm 3 Finding the real time counterterms and propagator simultaneously: 

1. Start with T.P{t;p) =0. 

2 . Loop over the three-momentum p: 

(a) Initialise p{t~0;p)~0, p{t — at] p) = at . 

(b) Solve the time-explicit ODE (I21D for p{t;p) with the thermal mass squared on the left 
hand side . [79] 

(c) Deduce F{t;p) using Fourier transforms and Eq. (I22D . [Isoll 

(d) Store p{t;p) and F{t;p) for the relevant t-range . 

3. Update the self -energy T,P{t;p) and adjust SM'^ to 

/>OC 

6M^ = - Y.P{t;p^Q). (CI) 



4. Iterate from step [2]. 

5. Repeat this procedure at an other temperature and with the corresponding thermal mass and 
use Eq. g38D to obtain 5\ and 6m?. 

Algorithm 4 Generating a far- from- equilibrium self- consistently correlated quantum state (modes are populated ac- 
cording to a given mode temperature T{p)): 

1. Start with a Y,P(t;p} as obtained in Algorithm (3). 

2 . Loop over the three-momentum p: 

(a) Initialise p(t^O;p), p{t ^ at] p) = at . 

(b) Solve the time-explicit ODE ( [2TD for p{t]p). 

(c) Deduce F{t]p) using Fourier transforms and 

i^(c.; p) = (l/2 + 1/ [e-/^(^^) - l] ) p) . (C2) 

(d) Store p{t,p) and F{t,p) for the relevant t-range . 

3. Calculate the self -energy and iterate from step [2]. 

Algorithm 5 Solving the nonequilibrium evolution starting from a self- consistent initial condition. 

1. Generate a translation invariant Fs{t]p) and ps{t]p) using Algorithm (|4P . 

2. Set F{tx,ty]p) = Fs{tx - ty]p) and p{tx,ty]p) = ps{tx - ty]p) for -imom < tx,ty < 0. Set the 
evolution time to t = 0. 

3. Calculate the self energies TiF,p{t,ty,p) for t — tmom < ij/ < ^ • 

4. Use the discretised Eqs . ( I51D - ( [52l ) to obtain F{t + 5t,ty) and p{t + St,ty) for imom < ^j/ < 

5. Use the newly calculated propagator values in Eq. ( I51D to obtain F{t + St,t + St) . 

6. Increase t by dt and repeat from step [3]. 
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APPENDIX D: DISCRETISATION OF THE ISOTROPIC PROPAGATOR 



It is very natural to start from a discretised action before defining the 2PI effective action. The straightforward 
lattice discretisation of a homogeneous system will then lead to a propagator G(ii, x) where a; is a lattice 3- vector. 
The isotropy in x reduces to a permutation symmetry between the components of x. Using also the reflection symmetry 
one can conclude that if the action is defined on a lattice with an even linear size N, the total number of variables 
that describes a propagator at a given pair of time coordinates is {N + 6){N + 4:){N + 2)/48. This discretisation has 
been used by several groups, starting with Ref. [i^ , and also in Ref . [63| ■ 

We now introduce a different kind of discretisation, where the rotation symmetry of the propagator is exact, and 
the number of independent variables in the propagator is N. This corresponds to a discretisation on the level of the 
2PI effective action, too. 

In many cases we have to calculate a loop integral, which is the convolution of two symmetric functions 



{f *9){p) 



fi\k\)gi\p^k\)^- 
(2Try p 



^kfik) 
Ztt 



,2^'^-^^'^^- 

p-k\ ^~ 



(Dl) 



(Here and in the following we use plain letters for the modulus of a lattice 3- vector.) This convolution is not defined in 
the strict sense without specifying a UV regulator. By fixing a momentum-space cut-off one also introduces Umklapp 
processes: if an outgoing momentum after a scattering process lies off the Brillouin zone, one maps this back to one 
of the modes inside. We motivate our cut-off prescription with numerical convenience. We would like to calculate the 
convolution in Eq. (jDl[) as a product co-ordinate space: 



if*9){x) = f{x)g{x), 



where the f{x) is the inverse Fourier transform of /(fc) in terms of 

dk 



xf{x) 







27r2 



kf{k) sinkx . 



(D2) 



(D3) 



We will have to discretise this inverse sine transform in momentum and co-ordinate. To avoid problems with the 
potential singularities at a; — > and — > we choose a discretisation without these points. For the infrared behaviour 
will assume that kf{k) — > as fc — > 0, which is certainly true in a massive theory. For algorithmic convenience we 
extend the momenum space functions beyond the cut-off A by the equation kf{k) — (2A — fc)/(2A — A:), this implicitely 
determines the Umklapp behaviour. 

The discretised sine and inverse sine transforms matching these boundary conditions are known as DST-II and 
DST-III formulae, respectively, these are available in many numerical libraries ^69j . We introduce a lattice spacing a, 
which is related to the the highest stored momentum by A = 7r/a, the space points where we will have to use Eq. (|D2p 
are Xn = a {n + ^) , with n = 0..N — 1. In momentum space we define the momentum grid as kj = (j + 1)A/N, with 
j = 0..iV — 1. Since on a one-dimensional lattice of physical size L the momentum space lattice spacing is 27r/L, the 
corresponding "box size" in our discretisation is i = 2aN. With these notations the transformation rules read: 



fk = 



N 
k+1 

1 



J2 ("+^)^"^^"(^ (n+l]{k + l] 

n—0 



2 

N-2 



N-1 



(D4) 



^(fc + l)/,sin(^ L + ij(fc+l) 



k=0 



(D5) 



As we have already mentioned, the zero momentum and the equal time propagators are not stored at all, so we 
need separate formulae for the volume integrals: 



d^k 
(2^ 



2(a7V)3 

N-l 



N-2 



k=Q 



/ d^xf{x)^ina^y2 U+- /„ 
•^0 n=o V ^/ 



(D6) 
(D7) 
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We used this kind of cut-ofF discretisation in our presented numerics. In the production runs the time-hke grid was 
typically four times finer than the spatial lattice. 
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